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Abstract 

In this article we show that separation of variables for second-order superintegrable 
systems in two-dimensional Euclidean space generates both exactly solvable (ES) and 
quasi-exactly solvable (QES) problems in quantum mechanics. In this article we propose 
the another definition of ES and QES. The quantum mechanical problem is called ES 
if the solution of Schroedinger equation, can be expressed in terms of hyper geometrical 
functions mFn and is QES if the Schroedinger equation admit polynomial solutions with 
the coefficients satisfying the three-term or more higher order of recurrence relations. 
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1 Introduction 



It is well known that A^-dimensional nonrelativistic quantum systems described by the Hamil- 
tonian 

--Y,^-^ +V{Xi,X2,...,Xn) (1) 

are integrable if there exist linearly independent and global integrals of motion I^, £ — 
0, 1, ..N — 1 and Tq — 7i, commuting with the Hamiltonian (1) and with each other 

[le,n]^0, [T^,X,] = 0, i,j^l,2,...N-l. (2) 

This particular class of integrable systems is called superintcgrable (this term was introduced 
first time by S. Ranch- Wo jciechowski in [1]) if it is integrable and, in addition to this, possesses 
more integrals of motion than degrees of freedom. The additional integrals Ck, commute with 
Hamiltonian 

[A,^]=0, A; = 1,2, ...TV, (3) 

but not necessarily with each other. If the number D of linearly independent integrals takes 
the value D = 2N — 1 (N the number of degrees of freedom) then the system is called max- 
imally superintcgrable [2] and it is called minimally superintcgrable if it has D = 2N — 2 
integrals of motions. Three examples of this kind have been well-known a long time, namely 
the Kepler-Coulomb problem, the isotropic harmonic oscillator, and the nonisotropic oscillator 
with commensurable frequencies. 

The existence of additional integrals of motion for these systems leads to many interesting 
properties unlike standard integrable systems. In particular, in quantum mechanics, there 
is phenomenon of accidental degeneracy when all energy eigenvalues are multiply degenerate. 
This property is intimately related to the existence of a dynamical symmetry group, a so-called 
hidden symmetry group, which contains properly the geometrical symmetry group describing 
this system. For instance, the dynamical group for the hydrogen atom is 0(4) for the discrete 
spectrum [3], and the Lorentz group 0(3,1) for the continuous spectrum [4]. For the isotropic 
harmonic oscillator it is SU(3) [5]. 

In classical mechanics the additional integrals of motion have the consequence that in the 
case of superintcgrable systems in two dimensions and maximally superintcgrable systems in 
three dimensions all finite trajectories are found to be periodic; in the case of minimally super- 
intcgrable systems in three dimensions all finite trajectories are found to be quasi-periodic 
[6]. 

One of the most important properties for superintcgrable systems is multiseparability, the 
separation of variables for the Hamilton- Jacobi and Schrodinger equations in more than one 
orthogonal coordinate system [7, 8, 9]. For instance, the isotropic harmonic oscillator in three di- 
mensions is separable in eight coordinate systems, namely in Cartesian, spherical, circular polar, 
circular elliptic, conical, oblate spheroidal, prolate spheroidal, and ellipsoidal coordinates. The 
Kepler-Coulomb potential is separable in four coordinate systems, namely in conical, spherical 
parabolic, and prolate spheroidal coordinates. 
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A systematic search for such systems in two- and three-dimensional Euchdean space was 
started in the pioneering work of Smorodinsky and Winternitz with coUaborators in [10, 11, 12] 
and was continued in [2]. Particularly, in [11] it was shown that in two-dimensional Euclidean 
space there exist four superintegrable potentials (see Table 1.), three of them could be considered 
as the singular generalization of Kepler-Coulomb, circular oscillator and anisotropic oscillator 
systems. These results were extended for two- and thee-dimensional spaces with constant 
curvature (both positive and negative) [13], and on the complex two-dimensional sphere and 
Euclidean space [14, 15, 16, 17, 18]. The same program is continuing nowadays both in spaces 
with constant and nonconstant curvature [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. 

The wide applicability of superintegrable systems, both in physics and mathematical physics 
has stimulated further investigations. In the last fifteen years superintegrable systems (in spaces 
of constant curvature including flat Euclidean space) became a subject of investigation from 
many point of view: in [13, 18, 31, 32] via the path integral approach, in [19, 21, 22] by solving 
the Schrodinger equation with the help of the Niven ansatz [33], in [34, 35, 36, 37, 38, 39] 
from the purely algebraic approach. As has been shown by a number of authors, almost 
all superintegrable systems (excluding the pure A^-dimensional Kepler- Coulomb and harmonic 
oscillator potentials) generate an algebraic structure which may be considered as a nonlinear 
extension of the Lie algebra (in classical mechanics Poisson algebras), namely a quadratic 
algebra. The general form of quadratic algebras, which are encountered in the case of two- 
dimensional quantum superintegrable systems has been investigated by Daskaloyannis [39]. 

In spite of all the above listed characteristics, superintegrable systems would not be so 
useful, except for the property of exact solvability of the all known superintegrable systems 
(at least in case of the second kind of supcrintcgrability). More precisely this means that after 
any separation of variables each of the separated ordinary differential equations admits an exact 
solution. This question, again quite recently has been discussed in the literature [40, 41]. As 
mentioned in these papers the term exact solvability is defined quite differently by different 
authors. ^ Indeed, in spite of an "intuitive" understanding of the term exactly solvable, no 
universal definition exists up to now. 

On the other hand there are limiting cases of well-known one- dimensional exactly solvable 
systems, namely the harmonic oscillator and Coulomb problems with (7 > —1/4) in- 

teraction, Morse potential, trigonometric and modified Poschl- Teller potentials, trigonometric 
and hyperbolic Manning- Rosen potentials [43, 44], and the Natanson potential [45]. All these 
potentials have the general property that the Schrodinger spectral problem for bound state (or 
continuous state) has an explicit formula for the whole energy spectrum, and the eigenfunctions 
(up to the asymptotic ansatz or gauge transformations [40, 41]) are of hypergeometric type 
iFi, 2F1. For the bound states we have solutions in term of classical polynomials [46] whereas 
for continuous states just infinite series. Moreover, hypergeometric functions describe both the 
continuous quantum systems as well as the finite systems and appear also as solutions of related 
difference equations, for instance, the finite one- and two-dimensional oscillator expressed in 
terms of discrete variables polynomials: Krawchuk, Meixner and Hahn [47] . 

Thus, we propose another definition of exact-solvability: a quantum mechanical system 

^In [40, 41] (see also the recent paper [42]) we read that "an exactly solvable quantum mechanical system 

can be characterized by the fact that in its sohition space one can indicate explicitly an infinite flag of functional 
linear spaces, which is preserved by the Hamiltonian" or the "Hamiltonian is exactly solvable if its spectrum 
can be calculated algebraically" . 
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is called exactly solvable if the solutions of Schrddinger equation, can be expressed in terms 
of hypergeom,etric functions mFn- (Basically, we are requiring that the coefficients in power 
series expansions of the solutions satisfy two-term recurrence relations, rather than recurrence 
relations of higher order.) It is obvious, that an A^-dimensional Schrodinger equation is exactly- 
solvable if it is separable in some coordinate system and each of the separated equations is 
exactly-solvable. Further, We say that a superintegrable system is exactly solvable if is exactly 
solvable in at least one system of coordinates. 

At first sight, such a definition of exactly-solvable problems may seem too narrow, but it 
leads us to distinguish two kind of models: 1) those which is possible to study analytically and 
2) those which admit just numerical solutions (even if they admit polynomial solutions and 
their energy spectra can be calculated by solving an algebraic equation or system of equations). 

The process of separation of variables in the N- dimensional Schrodinger equation leads 
to ordinary differential equations having as solutions many special functions of mathematical 
physics. A further complication of the separated equations involves the A^" separation constants. 
In general we have a multiparameter eigenvalue problem [48] . It is possible to distinguish three 
different cases, namely when there is complete, partial or non-separability of the separation 
constants. It is obviously that in the case of complete separability (of separation constants) the 
initial A^-dimensional Schrodinger equation splits into A^ independent second order differential 
equations, each involving a single separation parameter. This situation occurs, for instance, in 
the case of separation of variables in the Helmholtz (free Schrodinger equation, which is also su- 
perintegrable) or the Schrodinger equation for the harmonic oscillator in Cartesian coordinates. 
The second "extremal" case, when complete non-separability exits, is realized, in separation 
of variables for the same problems but in ellipsoidal coordinates. In the last case the each 
separated second-order differential equation contains at once all separation constants (usually 
depending from dimensional or non-dimensional parameters) [8, 9], for which the simultaneous 
quantization becomes nontrivial. 

The standard method of solution of a second order ordinary differential equation, obtained 
after separation of variables in - dimensional Schrodinger equations, involves (after taking 
into account the asymptotic ansatz) expansions around one of the singular points of the differ- 
ential equation (standard power-series method [49], or the so-called Hill-determinant method 
[50]). After that the problem reduces to solution of the recurrence relations for the expansion 
coefficients. If the coefficients obey a two-term recurrence relation, then the corresponding 
solution will be written in closed or analytic form or in terms of hypergeometric functions and 
we have an exactly-solvable problem. Such situations occur when separation of variables 
for superintegrable systems is possible in sub-group type coordinate (spherical, cylindrical and 
Cartesian) [51] and often in parabolic type coordinates. This method also powerfull when sepa- 
ration of variables possible in non-subgroup systems of coordinates as spheroidal or elliptic. In 
this case we arrive to high-order recurrence relation, the subsequent analysis of which, allow to 
investigate the behavior of the solution and to answer on the question whether exist the finite 
solution. 

Actually there is another general approach for solving the Schrodinger equation by exploring 
the Niven - type ansatz [33] and essentially is based on existing of finite solutions. According 
to this methods the complete solution can be constrated without direct separation of variables 
and computed in terms of zeros of corresponding polynomial. This method have been using in 
papers [19, 21, 22] at the investigation of two- and three-dimensional superintegrable systems 



3 



in Euclidean and curve spaces. 

Consider now the problem of motion in the plane for a charged particle with two fixed 
Coulomb centers with coordinates (±D/2,0) (so-called plane two center problem) 



0L2 



y2 + (x + D/2)2 ^y^ + {x-D/2f 



(4) 



This system is not superintegrable and separation of variables is possible only in two-dimensional 
elliptic coordinates (see eq. (92)). Upon the substitution ip{u, ji; D"^) = X^u; D'^)Y{fi; D'^) and 
the separation constant A{D), the Schrodinger equation splits into a system of two ordinary 
differential equations 



+ 



D^E 



D^E 



cosh^ u + D{ai + cosh u + A{D) 



cos^ fi + D{ai — ^2) cos fi + A{D) 



Y = 0. 



(5) 



(6) 



Both equations (5)- (6) belong to the class of non-exactly solvable problem. In generally the 
polynomial solutions do not exist even for the case of discrete spectrum < 0^, and each of 
the wave functions X{i';D^) and Y{iJ,;D^) is expressed as an infinite series with a three-term 
recurrence relation. 

Let us now put 0:2 ~ 0. Then the potential (4) transforms to the ordinary two-dimensional 
(2d) hydrogen atom problem, which is well-known as a superintegrable system [52, 53, 54] 
with dynamical symmetry group SO (3), and admits separation of variables in three systems of 
coordinates: polar, parabolic and elliptic. In this case we can see that the separation equations 
(5) and (6), namely 



fx_ 

fY 



+ 



D'^E 



D^E 



cosh^ u + Dai cosh u + A{D) 



X = 0, 



cos^ fj, + Dai cos fj, + A{D) 



Y = 0. 



(7) 



(8) 



transform into each other by the change ^ ii/. Thus separation of variables in elliptic coor- 
dinates for the 2d hydrogen atom gives two functionally identical one-dimensional Schrodinger 
type equations with two parameters: coupling constant E and energy A{D) (correspondingly 
energy and separation constant for 2d), but one defined on the real and the other on the imagi- 
nary axis. In other word, instead the systems of differential equations (7)-(8), the task reduces 
to solving only to the one of equations (7) or (8) for what the "domain of difinition" is complex 
plane. The requirement of finiteness for the wave functions in the complex plane permits only 
polynomial solutions (see for details [55]). As result we obtain simultaneously quantization of 
the energy spectrum 



En — 



2(n+ 1/2)2' 



n 



0,1,2,... 



(9) 



^To be completely correct let us note that the polynomil solution exist only for special values of parameters 

ai, OLi and R. 
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and the elliptic separation constant As[D) where s = 0,1,2, ...n (as a solution of an nth-degree 
algebraic equation). The polynomial solution difine by the help of finite series with the three- 
term recurrent relations for coefficients. They cannot be considered as exactly-solvable and 
maybe investigated only numerically. A similar situation occurs, for instance, in the case of 
the two-center problem in three-dimensional Euclidean space (the so-called prolate spheroidal 
radial and angular Coulomb wave functions) [56] and three-dimensional sphere (Heun wave 
functions) [57], where after eliminating one of the Coulomb centers the problems have reduced 
to superintegrable systems admitting only polynomial solutions. The presented above (and 
many others) examples let us to claim about deep connection of the notion of superintegrability 
and existing of polynomial solutions of the corresponding Schrodinger equation. 

At the other hand side the each of equations (7) or (8) have the form of one dimensional 
Schrodinger equation with the parameter E and for eigenvalue A{D), and could be separately 
considered in the regions e [0, 27i] or u e [0, oo) correspondingly. Then for the arbitrary 
values of constant E ^ the solutions of eqs. (7) or (8) expressed via infinite series and only on 
the "energy surface" corresponding of 2d hydrogen atom (9) the solutions split onto polynomial 
and nonpolynomial sectors (each of these sectors are non complete) and for fixed number n, 
only part of eigenvalue As{D), (s = 0, 1, 2...n) are possible to calculate as nth-degree algebraic 
equation. We can say that the eqs. (7) and (8) "remember" for its polynomial solutions. It 
is obviously that the spectrum of As{D), (s = 0, l,2...n) and solutions in polynomial form of 
the each of equations (7)-(8) coinside with the eigenvalue of separation constant and the wave 
function after the reduction to one of the regions n e [0, 27r] or u e [0, oo)] for 2d hydrogen 
atom. 

This phenomena have been intensively discussed in literature in the late 1980's and was de- 
termined as quasi- exactly solvability (this term first time introduced by Turbiner and Ushveridze 
in [58]) and such models called as quasi-exactly solvable systems [59, 60, 61] (see also [62] and 
references therein). The crucial example stimulated the investigation of quasi-exactly solvable 
systems is the hamiltonian (1) with the anharmonic potential 

V{x) = -u^x' + 2/3a;V + {2l3^u:^ - 2Su; - X)x^ + 2 ^^^"^^^ (10) 
2 x^ 

where u, f3, 5 > 1/2 and A are the constants. As it was notice of many authors [63, 64, 65], 
this systems only for special values of constant A = uj{2n + 1), (n = 0, 1,2...) admit partial 
polynomial solutions 

There are two different approaches in investigation of quasi-exactly solvable systems. In the 
base of algebraic approach formulated by Turbiner in [59] the fact of quasi-exactly solvabilily 
was explain in term of "hidden symmetry algebra" sl{2, R) ^. More precisely it means following: 
the one-dimensional Hamiltonian (1) after suitable changes of variable z — ^{x) and "gauge 
transformation" H — e~"^^^Ti.e"^^^ can be written in form 

CabJaJb+ E CaJa (12) 

a,6=0,± a=0,± 

^For example when = {n ^ oo) the equations (7) and (8) transforms to periodic and modifying Matie 
equations, which are non-exactly solvable. 

^Really this is not hidden dynamical symmetry in usual sence because the hamiltonian (12) belong to the 
enveloping algebra but not a Casimir. 
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where the first-order differential operators {J±, Jq} span to the commutation relations for al- 
gebra sl{2,R) [59]. 

The mentioned above analysis for the 2d hydrogen atom have shown, that in spite of the 
elegance of algebraic approach, the phenomena of quasi-exactly solvability have really more 
deeper roots that it can be explain via "one-dimensional" model (12). The other examples are 
hydrogen atom and oscillator problems on two and three-dimensional spheres [19, 66] and two- 
dimensional hyperboloids [22], which generate not only hyperbolic and trigometric but elliptic 
type of quasi-exactly solvable systems (se also oldest articles [67, 55, 68, 69]). We can also 
mention Lam'e polynomials. They comes from separation of variables for the Helmholtz (also 
superintegrable!) or Schrodinger equation in elliptic coordinates on the two-dimensional sphere. 
As it have been also determined in [37] (without showing of the mechanism of this phenomena) 
some of the quasi-exactly solvable systems can be obtained through dimensional reduction from 
the two and three-dimensional superintegrable models with quadratic invariants (second-order 
superintegrability) . 

The second approach, known as analitic, was formulated by Uschveridze (see for example 
[60, 61, 62]) and really represent a one-dimensional reduction of the Niven-Stilties method for 
solving of many-parametric differential equation as generalized Lame equation (or ellipsoidal 
equation) [33] . The solution in this method are presented in the term of zeros of polynomials 
Pn{x'^)- Then the wave function (11) can be rewritten in the form 

*(x) ^ x^^-'^ e-^^'-^'"^' n^=o(^' - Ci), (13) 

where the numbers (Ci, C2, •••Cn) satisfying the systems of n algebraic equations (see section 
2.3). Accordingly to the oscillation theorem the number of zeros in physical interval e [0, 00) 
enumerate the ground state and first n - excitations, describing in terms of all zeros (complete 
solutions of the systems of algebraic equations and including non physical section e (—00, 0]) 
as 



n 1 

1=1 



(14) 



The two nutural questions appering in such approach: what is the physical meaning of the 
negative zeros ^j, and why in the correct formula of energy spectrum (14) are participating all 
n zeros of the polynomial Pn{x'^)7 

With this article we begin the new investigation of all known two- and three-dimensional 
second order superintegrable systems on curved spaces (Euclidean, sphere, hyperboliod and 
pseudo-euclidean) base on the separation of variables and direct solutions of the Schrodinger 
equation. We pay special attention to non-subgroup type coordinates and prove the existence 
the polynomial solution for this systems, which we found surprisingly that never have been 
done before. 

We also demonstrate that the fact of quasi-exactly solvability are directly related with the 
multiseparability of the second-order superintegrable systems, from one hand side, and with 
presence of polynomial solution in this systems on the other hand side. 

The present work devoted only two (singular anisotropic and singular circular oscillators) 
from four famous superintegrable systems in two-dimensional Euclidean space (see first two 
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potentials in Table 1.). The next two systems maybe transformed (only for the discrete spec- 
trum!) to singular circular oscillator (for V3) or ordinary shifted oscillator (for V4) systems by 
the help of Levi-Civita mapping [70] . 



Table 1: Superintegrale potentials in two-dimensional Euclidean space. 



Potential V{x, y) 


Coordinate 
System 


1 1 1 
V,^'y{Ax' + y') + k,x+'''' 

2 2 


Cartesian 
Parabolic 




Cartesian 

Polar 
Elhptic 


\/x^ + y^ A\/x^ + y'^\\/x^ + y'^ + X \/x^ + y'^ — x J 


Polar 

Elhptic II 
Parabolic 




Mutually 
Parabolic 
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2 Singular anisotropic oscillator 



Let us first consider tiie potential (/c2 > 0) 

1 



^2 _ 1 

Vi{x, y) = ^u;\4x^ + y^) + hx + ^— ^ 
2 



(15) 



(second potential on the table 1.), which we will call the singular anisotropic oscillator. The 
Schrodinger equation has the form 



92 

+ 



dx^ dy'^ 



2E - uj^{4x^ + y^) - 2kix - 



r 



^ = 0. 



(16) 



For k2 > 1/2 the singular term at y = is repulsive and the motion takes place only on one of 
the half planes {—00 < a; < 00, |/ > 0) or {—00 < x < 00, y < 0), whereas for < /c2 < 1/2 in 
whole plane (x, y). There are two coordinate systems of relevance here: Cartesian and parabolic 
coordinates. 



2.1 Cartesian bases 

Separation of variables for the eq. (16) in Cartesian coordinates leads to the two independent 
one-dimensional Schrodinger equations 



dx'^ 



^ + (2Ai - Au^x"^ - 2kix)^i = 0. (17) 
+ (2A2-u;V-^MlV'2 = 0. (18) 



-— + [2X2 - u; y 

dy^ V y 



where 

*(x, y; fci, ±k2) = '0i(a:; A;i)^/'2(?/; ±k2) (19) 

and Ai, A2 arc two Cartesian separation constant with Ai + A2 = E. 

The equation (18) represents the well-known linear singular oscillator system studied in 
detail (see for instance the books [71, 72] and articles [11, 20, 73]). It is an exactly-solvable 
problem and has been used in many applications, for example as a model in N - body problems 
[74], or fractional statistics and anyons [75, 76]. The complete set of orthonormalizcd cigen- 
functions, (on 1/2) in the interval < y < 00 of eq. (18), can be express in terms of finite 
confluent hypergeometric series or Laguerre polynomials 



ipn2{y;±k2 



and A2 = a;(2n2 H- 1 ± ^2). We assume that the positive sign at the k2 has to be taken if 
k2 > \ and both the positive and the negative sign must be taken if < /i;2 < |, so that the 
polynomials have finite norm. Let us also note that unlike the potential (15) the wave function 
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is not invariant under the replacement /c2 —k2 and splits into two families of solutions that 
transform to one another under this change. 

The second equation (17) easily transforms to the ordinary one-dimensional oscillator prob- 
lem. In terms of Hermite polynomials the orthonormal solution (in region — oo < x < oo) 
is 

V'ni(:r;A;i)= (^— j -j===H^,{V2^z), z^x+^, (21) 

and where Ai = uj{2ni + 1) — g^- Thus the complete energy spectrum is 

E ^ Xi + X2^u;[2n + 2±k2\ - n = m + ria = 0, 1, 2, ... (22) 

and the degree of degeneracy for fixed principal quantum number n is (n -|- 1). Finally note 
that the separation of variables in Cartesian coordinates leads to two exactly solvable one- 
dimensional Schrodinger equations and the complete wave function may be constructed with 
the help of formulas (20), (21) and (19). 

2.2 Parabolic bases 
1.2.1. Sepeiration of variables. 

Parabolic coordinates ^ and rj are connected with the Cartesian x and y by 

x=\{e-v') , y = ^v , ^eK,rj>0. (23) 
The Laplacian and two-dimensional volume element are given by 

The Schrodinger equation in parabohc coordinates (23) is 



1 /(92^ <92^\ 



i-2 _ 1 

2^C4 _ t2^2 _i_ „4\ _ u (c'i _ _ '^2 4 

2^2 



2E-u\e-eri'-rv')-k,{e-ri') 



* = 0. (25) 



Upon substituting 



and introducing the parabolic separation constant A, the equation (25) split into two ordinary 
differential equations: 

^ + (2Ee-u^r-k^e-^^X = -XX, (26) 
+ (2Er^'-u;'r]' + k^r,'-'^^)Y = +XY. (27) 



9 



Equations (26) and (27) are transformed into one another by change ^ < — > irj. We have 

*(e, V, E, A) = C{E, A)Z(e; E, \)Z{ir^- E, A) (28) 
where C{E^ A) is the normahzation constant determined by the condition 

/ dr] d^{e + v')\'^{^,V;E,X)\'^l (29) 

Jo J-oo 

and the function Z{ii; E, A) is a solution of the equation 



Z{n;E,X) = XZ{n;E,X). 



(30) 



Thus, at ^ e (—00,00) we have eq. (26) and at e [0,ioo) - the eq. (27). Note that in the 
complex /X domain the "physical" region is just the two lines Im /i — and Re // = 0, Im > 0. 
Our task is to find the solutions of eq. (30) that are regular and decreasing as /i — > ±00 and 
II — > ioo . 

1.2.2. Recurrence relations. 

Consider now the equation (30). To solve it we make the substitution 



Ziii; E, A) = exp (^-^/ - ^^^j V(/^; E, A), 



and obtain the differential equation 



+ 



— 2ujijb jj, + 



2a;2 



dip 
d/i 



where 



2Ei/ + A ^ = 



ki 



E^ E+-^-u;{2±k2), X^X-—{l±k2). 



6uj 



Passing to a new variable z = /i^ in eq. (32), we have 



[1 ± k2) - U}Z [ z + 



2u;2 



dijj ^ 
dz 



1 It 
-Ez + -A 

2 4 



^ = 



We express the wave function ^l){z) in the form 



(31) 



(32) 



(33) 



(34) 



i;{z-E,X)^Y.ME,X)-^ 



(35) 



The substitution (35) in cq. (34) leads to the following three-term recurrence relation for the 
expansion coefficients Ag = As{E, A), 



(s + l)(s + l±fc2)A,+i + 



+ 



A - — (2s + l±A;2) 



E 



— a; (2s ± /C2) 



As-i = 0, 



(36) 
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with the initial conditions A_i — and Aq — 1. 
1.2.3. Asymptotic behavior. 

To understand the behavior of the coefficients for large s we use continued fractions theory [77] . 
Setting 



As 
As-i 

fis) 



U{s) (37) 

oj rd + i)r(±f + 1 + i)r(^ ± f + 1 + 1) 
2r(f + i)r(±f + f + i)r(^±f + | + 



where T{z) is the Gamma function, we can write the recurrence relation (36) in the standard 
form 

= (38) 



2 r(| + i)r(±f + f + i)r(|- ± f + f + i)[-a + ^(2. + 1 ± k,)] 



r(| + i)r(±f + 1 + i)r(f ± f + § + 1)23 

where a = -{E + ki/{8u;^))/2. Note that 



{S + I){±k2 + S + 1)' 



Stirling's formula for the Gamma function T{z) = ^/^e ^v^27r(l + 0(1/2;)) as \z\ oo with 
|arg z\ < 71, gives f{s) = -^f (1 + 0{l/s)) and bs = ±2;;^{^ + 0{l/s)). In the following we 
take ki < 0, /c2, A, £^ real and a; > 0. Without loss of generality we can assume bs is positive for 
sufficiently large s since, otherwise, we could make the replacements bs — > —bs, ~Cs- 

Since X^fcg = oo, it is a consequence of the Seidel-Stern theorem that the formal continued 
fraction expressions for the converge: 



Cs = bs+bs+i-\ \-bs+k+ 

Moreover, standard continued fraction theory tells us that 

where 



hm (39) 



and 



(40) 
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Furthermore the relation A^^b!^1i ~ ^i-i-^i*^ = ("1)" ^ holds for all n > 0, which implies 

Al!^ A«i (-1: 



^ 

d(s) r(«) ~ r(*') r(s) ' 



This result in turn implies that the sequence A2I/ is, for large s and n monotone increasing 
in n and goes to in the limit, whereas >l2n+i/-^2n+i monotone decreasing in n and goes to 
^5 in the limit. For example. 



4(5) a{s) 7 

^2n+2 ^2n 02^+2+3 



r(s) r(s) r(*) r(s) 

-"2n+2 J-'ln ^2n ^2n+2 



(41) 



> V-J^2n+2+^ . 



It follows from (39), (41) that 

£>0 n=l -D2n -"271+2 

Simple estimates using the recurrence relations (40) give 

n n 
^2n > 1 + X] ^2m+s, -620+1 > X] ^2m+l+s, 



m=0 



Substituting these results into the identities 



-"2n — (^2771+5-02771-1) -°27i+l — (^2771+3+1-02771 

777=1 771.=0 

we get refined upper bounds for -B2n > -S2n+2- We can approximate the sum YJ'm=s by the 

integral j'g ^ (ia; and use similar approximations to get an upper bound for the series (42): 

I . I dy 

Jo + s{y^ + 1) 

for positive constants Kj independent of s. This shows that is uniformly bounded in s. 
Since ^5+1 = —bg + a-nd 6^ — > as s — > 00 it is also true that |l/^s| is uniformly bounded 
in s. 

It follows from (38) that 

_ e.-i - (bs + ^s-i) (1 - bs-i^s-i) 

?s+l - ?s-l — 



1 - bs-iCs-1 

Now choose Sq so large that bg+i < bg and bg^s < 1 for all s > Sq. Note from this identity that if 
Csi-i > 1 for some Si > Sq then ^s^^+i > ^si-i > 1- Thus the sequence ^s-^+2fc-i is monotonically 
increasing for all k >0. Since \^s\ is bounded, it follows that in this case limfc^oo Csi+2fc-i = C+ 
exists, and ^+ > 1. Since ^^+1 = —bg + 1/^^; bg — > as s — > 00 and \l/is\ is uniformly bounded 
in s, then the sequence ^si+2k is also convergent, lim^^oo ^si+2fc = ^- where < ^_ < 1. 

The other possibility is that ,^ < 1 for all s > sq. Since 1/^^ — Cs+i = bg ^ as s —>■ 00, and 
> 1, C,s+i < 1 for all s > So it follows that limfe^oo^fe = C+ = C- — 1- Thus in all cases the 
sequences .^2*; and ^2A;+i converge. 
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We conclude that 



l + 0(l/s)), ^+^ = 1, 



depending on whether s is even or odd. Thus asymptoticaUy As ~ v^yTf) depending on 
whether s is even or odd, and 



(43) 



Then we have for z > [the case of eq. (26)] 




SI 



— ^± cosh ( 



+ sinh ( -z' 



(44) 



This function does not belong to the Hilbert space. If /ci > then we must make the replace- 
ments bs —bs and —Cs- This has the effect of replacing z by —z in (43). Now the 
asymptotic solution is oscillatory. However, then for 2; < [the case of eq. (27)] the solution 
doesn't belong to the Hilbert space. The solution we have found is the minimal solution of the 
thrcc-tcrm recurrence relations. There is a linearly independent solution, but the coefficients 
grow more rapidly than the minimal solution coefficients. 

1.2.4. Energy spectrum and sepEiration constant. 

Thus the function Z(iJ,) cannot converge simultaneously at large n for real and imaginary 
and therefore the series (35) should be truncated. The condition for series (35) to be truncated 
results in the energy spectrum giving the same formula (22) and now the coefficients Ag = 
74"^(A;i, ±k2) satisfy the following relation 



(s + l)(s + 1 ± k2)As+i + psAs +uj{n+l- s)A_i = 0, 

A ki 



(45) 



/3. = --^(2. + l±M. 

The three-term recurrence relation (45) represents a homogeneous system of n + 1 - algebraic 
equations for n -|- 1 - coefficients {Aq, Ai, A2, ...An}. The requirement for the existence of a 
non-trivial solution leads to a vanishing of the determinant 



Po l±k2 

un (3i 2(2 ± k2) 



2UJ Pn-l ^{71 ± k2) 



(46) 



The roots of the corresponding algebraic equation give us the (n+1) eigenvalues of the parabolic 
separation constant A„(/ci, ±A;2). It is known that all roots for a such determinants are real and 
distinct [78]. Thus all values of the separation constant are real and can be enumerated with 
the help of the integer q, namely the values are A„(A;i, ±^2) ^nq{ki, ±A;2), where < q < n. 
The degeneracy for the n - energy state, as in the Cartesian case, equals n + 1. 
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Note that eq. (46) is invariant under the simultaneous transformation ki — > —ki and 
A —A. Thus if one of the A = A„(A';i, ±A';2) is root of cq. (46), then A = — A„(— A:i, ±A';2) 
is also a root of the same equation. We see that for the odd energy state (n-odd) the range 
of Xnq{ki,±k2) sphts into two subset A^^^^ and A^^^ connected with each to other by the rela- 
tion A^^J(/ci, ±A;2) < — > — A|fJ,(— /ci, ±/c2)- For n-even, always there exists the additional root 
Xnq{ki, ±A;2) = —Xnq{—ki, ±A;2), which equals zero when ki — 0. 

1.2.5. Wave functions. 

We will term the polynomial solutions of eq.(34), or eq.(32), as Mknq{z; ki,±k2), and the 
function (31) as Tanq{z; ki, ±k2) ^. Then the physical admissible solutions of eq. (34) have the 
form 

n 

Mknq{z;k^,±k2)=i^{z;E,X) = Y,Al'ih,±k2)z\ (47) 

s=0 

and the corresponding solution of eq. (31) is 

Tanqii^; ki, ±k2) = exp (^—/^^ - ^/x^j /x^^'^^ Mk^qiiJ.^; ki, ±k2) (48) 

Observe that parabolic wave functions (as also Cartesian wave functions) split into two classes 
and transform to each other via k2 —k2. In the case A;2 = (when the centrifugal term 
disappears), the solution (48) becomes an even and odd parity wave function under exchange 

It is known that there exists a direct connection between the quantum numbers q and num- 
bers of zeros of the polynomial (47) and, therefore, the eigenvalues of the separation constant 
A„g(/ci, ±^2) may be ordered by the numbers of nodes of the wave function Ta„g(yu; fci, ±^2)- 
Indeed by we will see that these are orthogonal polynomials, hence [49] , all the n - zeros of the 
Mknq{z;ki,±k2) are situated on the real axis —00 < z < 00, and all zeros have multiplicity 
one. Assume that the separation constants Xnq{ki, ±/c2) are numerated in ascending order, i.e 

A„o(A:i, ±k2) < Xni{ki, i/ca) < < A„,„_i(A;i, ±k2) < Xn,n{ki, i/ca) (49) 

then according to the oscillation theorem [79] , the quantum number q also enumerates the zeros 
of polynomials Mknq{z; ki, ±^2) in the region 2; > 0, or the real axis of fi, (sec cqs. (52) - (60)). 
Let us now introduce two quantum number qi and q2, which determine the zeros of polynomials 
Mknq{z; ki, ±^2) for 2; > and z < 0, correspondingly. Then qi + q2 = n, and 

Xnqiiki, i/ca) = -Xnq^i-ki, i/ca) (50) 

For Ijl — ^ the function (48) gives the solution of equation (26), and for ji = ir] the solution for 
equation (27). Thus the parabolic wave function (28) can be written in following way 

*ngig2(?>^; ±^2) = C^q^q^iki, ±^2) Ta„gi(C; fci, ±^2) Tttnq^iW, h: ±k2). (51) 

1.2.7. The particular cases. 

^The notation Ta - devoted to the memory of Professor V.Ter-Antonyan (1942-2003) 
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Let us consider some low energy state n — 0,1,2. In case of n = we have 

Mkoo{z) = l, Aoo = — (liM (52) 
For n = 1 we have that g = 0, 1 and equation (46) admits two solutions 



Aio = -(2 ±k2)-J^ + 16u;(l ± k2). An = -(2 ± + J ^ + 16a;(l ± k2). (53) 
Therefore there are two polynomials 



h - \ kl + 16a;3(l ± fca) 
Mfcio(;^;A;i,±A;2) = 1 ^^^^^^^^j (54) 

Mku{z;k,,±k2) = 1 ^ ^^^^ ^ (55) 

For n — 2 the equation (46) is equivalent to a cubic algebraic equation 

A3 _ + [tt^ - 2u;(3 ± A;2)]A + 2A;i(l ± k2) = 0. (56) 

ZUJ lOJ^ 

At /ci = this equation can be simplified, and we have following three solutions 



A20 = -^32cu(3±fc2), A2i = 0, A22 = y3M3^M (57) 



with 

M*,(.) ^ l+\^gp.+ ^.= (58) 

M*.W ^ + (60) 

1.2.6. Orthogonality relations and normalization constant. 

The wave functions (51) as eigenfunctions of Hamiltonians are orthogonal for quantum number 
n, or for n' 

/ dr) dae + v')K',,J^,mh,±k2)K,,J^,mki,±k2)^o. (61) 

Because the energy spectrum is degenerate there exist additional orthogonality relations for 
quantum number q. Using the equations (26) and (27) it is easy to prove that for qi ^ q[ and 

?2 7^ ?2 

/CXD 
d^Tal. ii; ki, ±k2) Tan,,{i; h, ±/c2) = (62) 
-00 

roo 

/ driTa*^g,^{iri; ki, ±k2)Tanq2{w, h, ±^2) = 0. (63) 

J 
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Thus we have for 4 

/'■DO 



(64) 



Let us now calculate the normalization constant Cnqiq2{ki, ±k2). Prom the explicit form of the 
wave function 'i>^g^g^{^,'rj;ki,±k2) and the normalization condition (29), it follows that 



1 " 

-\Cr,,,,,{k„±k2)\' E i-lY'-'' A^{k^,±k2)A^,{k„±k2) 



s,s',t,t'=Q 



where 



2u; 



^±1/4 ^ i I 2 ^ 4 =t 4 



m=0 



m! 



(65) 



(66) 



2.3 Niven approach 

Let us express solution of the Schrodinger equation (16) in the following form [19] 



(67) 



From the eqs. (20), (21) and (31) follows that the function $(a;, y) is polynomial (product of 
two polynomials) from variables (x, y"^) in Cartesian coordiantes and (^^, rf') for parabolic ones. 
It satisfy the equation 

n^{x,y) ^ -2E^{x,y), 



where the operator is 



7^ 



92 92 
+ ^ + 



(1 ± 2k2) 



dx^ dy^ 
Taking into account that 



- 2a;|/ 



d_ 
dy 



X + 



A' 



(6^ 



(69) 



Mkng(z; ki, ±k2) = £ A^'^iki, ±k2) = n^=i(z - ae), 



(70) 



s=0 



where a^, i = 1,2, ...n are zeros of polynoms Mknq{z) on the real axis —00 < 2; < 00, and that 
in parabolic coordinates 



a a 



(71) 



we can choose a solution of eq. (68) in the form 



^x,y)^Mknq,{e;h,±k2)Mknq,{-v';ki,±k2)^U'l^, \y+2x-ae\ . (72) 
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Then from (68) follows that zeros ct^ must satisfy the systems of n - algebraic equations 

2 fl±A;2) h 



H uja£ 

on 



2a;' 



£= l,2,....n, 



(73) 



and for energy spectrum we again have a formala (22). The system algebraic equation (73) 

contains n - set of solutions (zeros) (a^^'*, ol2^ ^ OL^n)^ q — 1,2, ...n and all zeros are real. The 

positive zeros a£ > define the nodes of wave functions for equation (26), whereas negative 
zeros ae < define the nodes of wave functions for equation (27). 

The eigenvalues of parabolic separation constant can be calculated the same way via the 
operator equation A$(,t, y) = \^{x,y) (see for details [19]). More elegant way is direct to use 
the differential equation (34) [62]. Rewriting first the eq. (34) in following form 



' (P 



2c^2 



d_ 

dz 



A'mjz-—{l±k2) 



Mknq{z;kx,±k2) 



^XMkng{z;ki,±k2) (74) 



Putting now the wave function Mknq{z;ki,±k2) in form of (70), we arrive to the following 
result 



A„,(A;i,±A;2) =4(1±A;2) 



- + E 

4a; ^ 



1 



=1 a 



(9) 



(75) 



(in case of n = the sum must be eliminated) where the quantum number q — 1,2, ..n labeled 
the eigenvalue of parabolic separation constant. 



2.4 Interbasis expansions between Cartesian and parabolic bases 

We determine the interbasis expansion relating Cartesian and parabolic wave functions at the 
fixed energy 

ni+ri2 

^nq,q,{C,V;±k2)^ E (±^2) *n,n. (x, t/; ±A;2) . (76) 

ni=0 

For simplicity we consider only the case when ki = 0. To calculate the interbasis coefficients 
^nq^^ we use asymptotic methods [80]. First, we change in ^^g^gjlC; ^; ±^2) from parabolic 
coordinates to Cartesian ones 

= \J + y"^ + X, rf = \Jx'^ + y"^ — x, (77) 

and then let y tend to zero on the both sides of (76). As a result 

^ — > \x\ + X, ff' ^ — X, (78) 

and the dependence on variable y in (76) is removed. Using the orthogonality conditions for 
Hermite polynomials we get the equality 



W:^Zi^k2) = C^q,qA±k2) (^) ^ ^ ...j/'j^^l^^^^, 't,,) ,.(±^2) (79) 
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where 



/oo „ 
-00 

= r [Mknq^{2x; ±k2) + {-irMkr^g,{2ix; ±k2)] e"''"^' H^.iV^^x) dx. 

•J 

With the formula (47) the latter integral can be expressed via the coefficients A^{±k2): 

1 " r 1 / 2 \ * /"f^ 9 

^n^,..(±^2) = ^ E [1 + ^"(±^2) (-) z'^ H^,{z) dz. (80) 

It is more convenient to calculate L]!^^_^g^{±k2) independently for even and odd rii. By integration 
by parts we find that integral in the sum is nonzero only at ^ < s < ni + ^2, for ni - even and 
^^^^ < s < ni + ^2, for ni - odd. By changing the sum index and making the integration we 
finally get 



w:^Z^^\±k2) = c„,,,,(±^2) U:- 



TT \ 4 



T{n2 + I±k2] 



,n — 



2ujJ \ 2a;(i+"i±fe)(ni)!(n2)! 



X 



An r(2s+ni+2) 



for n even and odd respectively. 



(81) 



3 Singular circular oscillator 

The second potential of the singular oscillator is {ki, k2 > 0) 



V2ix,y)^lu;\x' + y') + l-'''' ^ ' 



-2 1 i.2 _ 1 ■ 
4 _^ '^2 4 



2" ■ " ^ ■ 2 V x^ 



(82) 



The corresponding Schrodinger equation separates in three different orthogonal coordinate sys- 
tems: Cartesian, polar and elliptical coordinates. 

3.1 Cartesian bases 

Let us first consider the separation of variables in Cartesian coordinates. From the asymptotic 
ansatz 

^(x, y) = x^^'^i y^^''^ ex.p[-u;{x^ + y'^)] X{x)X{y) (83) 
we obtain two independent and identical separation equations 

X{xi)^2XiX{xi), i^l,2 (84) 



+ I -2a; + ^ X- (1 ± 2ki)u; 



dz 



x: 



dxi 
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where Xi = x,X2 = y and Ai + A2 = —E. As in the case of singular anisotropic oscillator we 
assume that the positive sign at the ki has to be taken if /cj > | and both the positive and the 
negative sign must be taken if < /cj < |. 

The last equation is exactly the equation for confluent hypergeometric functions. The 
quantization rule gives 

\^ -uj{2ni±ki + l), Hi = 0,1,2,.... (85) 

and the solution of eq. (84) in terms of Laguerre polynomials X[xi) = L^j'^^uxJ). Thus the 
corresponding set of orthonormal eigenf unctions which are normalized in quadrant x > 0,y > 
(on 1/4) is 

*£,nf = c-Sv^^^) (x)^±'=Hy)^^'^e-t(-^+^^)L±^(^^')i^n^(^y') m 

where 



ni,n2 



(87) 



\\ r(ni ±ki + l)r(n2 ±k2 + l) 
Prom (85) we have 

En = uj{2n + 2±ki±k2), (88) 

where n = 711 + ^2 = 0,1, 2,... is principal quantum number and the degree of degeneracy is 
n+1. 

3.2 Polar bases 

The separation of variables in the Schrodinger equation for potential (83) in polar coordinates 

X = r cos 0, y = r sin 0, < r < 00, < (p < 27i (89) 
gives us the orthonormal solution in polynomial form 



"-"^ ^ ' ^' V r(n. + 2m±k^±k2 + 2) ' 

X e-'^L2f^i^'=2+^(cur2)$(^^i'±'=2)(0)^ n^, m = 0, 1, 2, ... (90) 



(2m ±ki±k2 + l)q\T{q ±ki±k2 + l) 



^| 2r(m ± A;2 + l)r(m ± fci + 1) 

X (cos0)V'±'=i(sin0)i/2±fep^±fc2,±fci)(^Qg20) (91) 

where P^'^\x) is a Jacobi polynomials and E = a;(2n ± /ci ± A;2 + 2), with n = + m and 
with the same degree of degeneracy (n + 1). 

Thus we have seen that quantum system (82) is exactly-solvable in the Cartesian and 
polar systems of coordinates. 
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3.3 Elliptic bases 



3.3.1. Separation of variables. 

The elliptic coordinate {u, jj) connected with Cartesian one by(0<i/<oo,0</x< 27r) 



D 



X 



cosh i/ cos//, y 



D 



sinh v sin 



(92) 



where D is the interfocal distance. As -D ^ and D oo, the elliptic coordinate degenarate 
into the polar and Cartesian coordinates 



cosh 1/ - 

sinh V cos /i 

The Laplacian and volume element are 



2r 



COS/i — > COS(f), 

2x 



D 



{D ^ 0), 
(D^ oo). 



+ 



(cosh 2u — cos 2/x) y du^ d/j,'^ j 
The Schrodinger equation with (82) can be rewritten as 



D2 

dV = — (cosh2i/ — cos 2 iji)di'djji. 
8 



(cosh 2v — cos 2^) —(cosh 2v — cos 2^) 



4 



COS^/i 



+ 



sin^ 



64 

(^1 ~ i) (^2 ~ \, 



sinh i/ cosh v 



= 0. 



and after the separation ansatz 

^P{u,|l■D^)=X{u■D^)Y{,,■D^) 



(93) 
(94) 

(95) 



(96) 



(97) 
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transforms to two ordinary differential equations 



+ 



4 



64 
64 



cosh 2i/ 



cos 2// 



cosli 2i/ 



^2 



+ 



sinh 1/ cosh v 



X = -X{D^)X, (98) 



O ^1 

cos^ 2// + 



cos^ H 



+ 



Y = +A(L>')y, 



(99) 



where A is the eUiptic separation constant. These equations can be written in the unit form 



+ 



64 



cos^ 2C 



D'^E 



■ cos 2^ 



1.2 _ 1 

4 



cos^ C sin ( 



Z{0 = KD')Z{C) (100) 



where at C £ [0, 27r] we have the equation (99) but at C £ [0, zoo) - equation (98). In other 
words, as we see from Fig. 1 , in the complex plane C, the "physical" are only the shaded domain 
on the two lines Im = and Re C = 0. 

For ki^2 > \ the centrifugal barrier is repulsive and motion takes place in only one of 
the quadrants, as ^ G [0,7r/2], whereas for < ki^2 < | it takes place in the whole region 
C e [0, 27r]. For the particular case ki — k2 — \ the equation (100) transforms to the problem 
of the ordinary two dimensional oscillator and has been investigated in detail in the paper [68]. 
In this article have shown that the solution of eq. (100) (for ki — k2 — 1/2) is, described by the 
Ince polynomials [81]. 

In the case where k\ and /c2 are integers, eqs. (98) and (99) coincide with those that have 
been found via separation of variables in the Schrodinger equation for the four dimensional 
isotropic oscillator in spheroidal coordinates [69] . 



3.3.2. Recurrence relations. 

Let us now consider the equation (100). First, introducing the function VF(C; -D^) according to 

cos 2C 



Z(C;L'^)=exp 
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we have the equation 



H — sm2C^^ + 



de 



dC 



D^u D^E 2, kf 
— - — cos 2C cos C — 



^2 



cos^ C sin C 



-A 



(101) 



W = 0. (102) 



For kl — k2 — 1/2 this is the Ince equation, [49]. 
Next the substitution 



W{C;D') = (sinC)^^'^^ (cosC)^^'^^ U{C;D') 



yields the equation 



+ [(1 ± 2A;2) cot C - (1 ± 2ki) tan C + ^ sin 2C] ^ + [p cos" ( - A]C/ = 0, 
dC 4 dC 



where 

p=^[uj{2±ki±k2)-El A = A + ^(l±fci) 



(1 ± A;i ± A;2)' - 



D'E 
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(103) 



(104) 



■. (105) 
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Passing to a new variable t — cos^ ( we find 

" + {^^ ^ ^'^^^ - ^) - (1 ± k2)t + ^t{t -^)]^+\\pt- m = 0, (106) 

Finally, looking for the solution of the last equation in the form 

oo 

U{t;D')= J2MD')f, (107) 

s=0 

for coefficients As{D^) we have the three-term recurrence relation 
(s + l)(s + l±A;i)^+i - 



s{s + l±ki± k2) + + J 



+ ^\p + D^uj{s-l)]A,_i^O, (108) 
with Ai = and initial condition Aq = 1. 
3.3.3. Energy spectrum and separation constant. 

In analogy with our asymptotic solution of the recurrence relation for the singular anisotropic 
operator in the parabolic basis we use continued fractions. For the minimal solution of the 
recurrence relations we find for <^ 1 

^^-^^ '''"^ + 0(^)V (109) 



Thus we have 



and 



l7(cosC) ~ E ^^cos^' C ~ e*P(^ <=os2C). (110) 

Therefore we see that for this case the function Z(cosC;-D^) as C ^ 'oo is not normalizable. 
There is a linearly independent solution of the recurrence relations, but the coefficients grow 
even faster. Hence it follows that the scries (107) should be truncated. The condition for series 
(107) be truncated give us already known formula for energy spectrum (88) and reduce to the 
polynomials: 

n 

s=0 

where now the coefficients As = A^^'^^'^^'^^^D'^) satisfy the following three-term recurrent rela- 
tions 

{s + l){s + l±ki)As+i + psAs — (n-s + l)A,_i = 0, s = 0,l,...n (112) 
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with 



1 

4 



1)2 

(2s + 1 ± /ci ± /ca)^ + ^^(2s -n + 4:±4:ki 



{2±k,± k2) - + \{D^) 



64 



(113) 



and A_i = ^n+i = 0. 

The recurrent relations (112) becomes a system of (n + 1) hnear homogeneous equations for 
the coefficients ^45. Equating the corresponding determinat to zero 



Dn{\) = 



/5o 



n 



Pi 



2{2±ki) 



n{n ± ki) 



= 



(114) 



leads to the algebraic equation of a (n + 1) degree which determine the eigenvalues of elliptic 

separation constant X^^'^'^'^^'^'^^D'^). The quantum number q — 0, 1,2, ..n enumerate {n + 1) 
roots of eq. (114) and therefore the degree of degeneracy as in the polar and Cartesian cases for 
n-th energy state equal n + 1. It is also known at the corresponding numeration the quantum 
number q define the numbers of zeros of polynomial (111), which has exactly n - zeros situated 
in the open interval < i < 00, and therefore, the elliptic separation constant X!;^''^''^''^\D'^) 
may be ordered also by the numbers of the nodes of the eigenfunction of equation (100). 

3.3.5. Wave functions. 

Thus the the condition of finitness of the solution of eq. (102) allows the following polynomials: 

n 

2<±fci,±fe)(^.^2) ^ (sinC)5±fe2 (cosC)^^'^^ ^Api'±'=2)(L>2)(^Qg^)2.^ ^^^5) 



while the corresponding solution of eq. (100) is 



nq 



(116) 



We will term the polynomials T^''^''^''^\(; D^) as the associated Ince polynomials. In the case of 
ki — k2 — 1/2 these polynomials transforms to the four type of the ordinary Ince polynomials, 
which are even or odd with respect to the changes C. —>■ — C and C — C + [81; 68]. 

At ( = fi the wave functions (116) give us the solution of angular equation (99), and for 
C, = iu the solution of radial equation (98). For the each of wave functions, radial or angular, 
corresponds definite numbers of zeros which could be presented by two quantum numbers qi 
and obeying the condition qi-\- q2 — n. Then the complete elliptic wave function (97) may 
be written as 



'^V,/^;^') = a,i,.(±fci,±fc2;i^^)^S^'^^^)(/.;i^^) z^^^^^^^^\iv- D% (117) 
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where Cngj^g^{±ki,±k2; D^) is the normahzation constant. It could be calculated from the 
condition 

r)2 roo 1 

— d,yj^ rf/x(cosh^^-cosV)^iS^'^^*(^'/^;^')*iS^^^^^ (118) 

In same cases it is more convenient to have alternative form for the wave functions (117). 
Equation (100) represents the eigenvalue problem for the separation constant A(D^) and the 
corresponding eigenfanction. It easy to see that the operator in the left-hand side of eq. (100) 
is invariant under the simulteneous transformations 

C^C+|, D'^-D^, h^k2 (119) 

and therefore the sabstitution — )• — and ki < — > k2 does not change the set of eigenvalues 
of the separation constant X{D^), but change only their numeration: 

As a result transformation (119) transform the solutions (116) into each other 

zS'''^'^\C;D') ^S^'^'^'(C + (121) 

The coefficients A^^'^^^^^^^-D^) are calculated from the three-term recurrent relations (112) 
after the changing D"^ — ^ —D^ and ki < — > k2. From (121) we conclude that elliptic basis 
(117) can be represented in the form 

'^n^, f^; D') = C„,,,,(±A;i, ±k2-. D') Zi%'^^^'^\f,; D') Zi%'-^'^^\iu + |; -D% (122) 

where the quantum number qi, q2 {qi + q2 = n) again have the meaning of zeros of the angu- 
lar and radial wave functions and the constant Cnqig2{^ki,±k2; D"^) satisfy the normalization 
condition (118). 

3.3.6. Orthogonality relations. 

The wave function (117) as eigenfunction of the Hamiltonians are orthogonal n ^ n' 

r lo ^')*Sfi^'^H^, D')dV = (123) 

Eqs. (98) and (99) allow to prove the property of double orthogonality for wave functions 
4f^'^'^nC;^'), namely 

r '^'^^*(^-; D^)Ztg^-^''^\iu; D^)dv = (124) 

J 

/' 4?"^'^^*(/^; ^^)^£^'^^^H/^; ^')rf/^ = o (125) 

J 

for q\ 7^ g'^ and 52 7^ ^^d therefore when q^ 

j^dv JJd^i {cosh' u- cos' ^i)^^t,tt''^^*i^^ (126) 

24 



3.3.7. Limit D^O. 

Now we will show that in the limit D ^ the elliptic bases transforms to polar one. At the 
limit D — > in eq. (114) it is possible to neglect the terms depending of D^, then 

DniX) = PoPi f3n = (127) 

Let one of the (5s &t s = m equal zero: (5m = 0. Then A(0) = —(2m + 1 ± fci ± ^2)^ and 
therefore at the limit D — > the quantum numbers qi and q2 transforms to m and n — m — rir 
correspondingly. For s ^ m, we have 

(5s{D^) — > (5s{0) = -{s-m){s + m+l±h± fcs) (128) 

It follows from (128) that the three term recurrent relation (112) at D — > split to two term 
recurrent relations 

{s + l){s + l±ki±k2)As+i+psAs = 0, s = 0,l,....m-l (129) 

and 

(5sAs ^(n - s + l)^^-! = 0, s = m + l,m + 2,...n (130) 

Prom two-term recurrent relations (129) and (129), we get 

A ( ^^2^ g^p {-m)s{m + 1 ± fci ± fca). , . 

^'^^ ^ s\{\±k,)s ' ^^^^^ 

for s = 0, 1, ....m — 1 and 

Am^siD ) (-1) ,^_^^^ + 2±k,±k2)s{l±k,U \ ' ^ ^ 

for s = 1, 2, ...n — m, and {m)s = r(m + s)/r(m). Then according to eq. (93) 

n 

Z(Jf-±'=^)(/x;D2) ^ (cos0)^±^Msm0)^^'=^E^^P')(cos0)'^ (133) 

s=Q 

^^....2) (,,.^2) ^0 (,)(!-.) (I) E^(2r)^^ (134) 

Substituting here the eqs. (131) and (132), we get 

Z^^^^'^^^\^Ji■ D^) ^ (cos(/))5±*^i (sin(/))5±'=2 ^Fi{-m, m + 1 ± ki ± k2;l ± ki, cos^ 0) 

(cos0)5±'=i (sin(/))^±^2 P^''^'^'''\(f)) (135) 



{-!)"' ml 



(l±fcl)r 

Zi%'-^'^^\iu;D') ^ (ifrn^±k,) {m+l±k,±k2)m 

(1 =t Kijm 
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X ( iFi{-nr;2m + 2± ki± k2;ujr'^) 

. .(2m+i±fe2) n^! (m + 1 ± fci ± A:2)^ 

{l±h)m{2m+l±ki±k2)m 
2 / 9„\ 2m+l±fei±fc2 

X e-^ (I) L^r^^'^^'^K') (136) 

Prom these formulas and normalization condition (118) we finally get that te eUiptic basis (117) 
at the limit D — > becomes the polar basis 

V, D') ^ ^S'^'^^V, </•) (137) 

3.3.8. Limit D oo. 

Let us now to investigate the Cartesian limit. At D — > (X) we can neglect all final terms and 
again the determinat (114) transforms to product of diagonal elements. Let now one of the /3s 
for s — rii equal zero: = 0, then 

ASf-^^^n^^) XS^'^''\D^) = ^ - ^i^n, -2n + 6±7k^T h) (138) 



here ni = 0, 1, 2, .... Analogically it is easy obtain 

AS-±'^^)(-D^) = ^ + ^(4n2 - 2n + 6 ± 7fc2 T ^i) (139) 



with n2 = 0, 1, 2, and from the condition (120) follows that rii + n2 = n and rii, ^2 are the 
Cartesian quantum number. From the eqs. (138) and (138), and (113) it is easy to get (s ^ rii, 
and s ^ 712) 

/3s{D') ^-=^ D^s - m), Ps{-D') ""-^^ -DMs - n2). (140) 

These formulas and the condition A_i = shows that in the limit D ^ 00 the three-term 
recurrent relation (112) transforms to two-term recurrent relations 

{s + l){s + l±ki)A,+i{D^) —{s-ni)A,{D^)=Q, < s < rii - 1 (141) 

{s + l){s + l±k2)As+i{-D^) + ^{s-n2)As{-D^) =0, < s < ^2 - 1 (142) 
Analogically using the condition A^+i — 0, we obtain 

{s-ni)A,{D^) + ^{n-s + l)A,^i{D^) = 0, ni + l<s<n (143) 

{s-n2)As{-D^) + ^{n-s + l)A,_i{-D^) = 0, n2 + I < s < n (144) 
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Prom the two-term recurrent relations (141) - (142) we have 

Am+s{D ) — > j , 1 < s < n - ni (146 

4^ s! 



and the same for Agl—D"^ 

)2\ D^oo 



2. D^oo Ar^^{-D'^) {n2-n). 



A^,U-D') ™ ' l<s<n-n, (148) 

We can investigate now the hmit D — > oo in the wave functions. Taking into account the 
conditions (94) and formulas (145)- (148), we obtain the following results 

4r f|)'"'i:^£^55^(4x^)' 



2x\ - 

-^j iFi{-ni-l±ki-ux'), (149) 

= (^) ii^i(-n2; 1 ± k2, uy''). (150) 
Using the connection between hypergeometrical functions and Lagerre polynomials [46] 



^n(^) = ^^^^iF,{-n- 1 + «; x). (151) 



finally we get 



X (2a;) (27/) ^±^=2 ^±fci {ujx^)L^^^ (cuy^) (152) 

which after taking the limit in the normalization constant coinside with the Cartesian basis of 
singular oscillator. 

4 Conclusions and summary 

Let us here to summarize the presented investigation of sup erintegr able potentials Vi and V2. 
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Wc have determined that solution of Schrodingcr equation for the potentials Vi maybe 
constractcd via separation of variables with two different ways. One of them is exploting the 
separation of variables in Cartesian coordinates, which lead two independent exactly-solvable 
equations (17) and (18), the each of them represent the one-dimensional non-parametric spec- 
tral problem where the cartesian separation constants Aj play the role of energy. It admit to 
get the solution in form of Lagerre and Hermite polynomials, quantize both separation con- 
stants and as result to obtain energy spectrum for two-dimensional Schrodinger equation. The 
second way is the separation of variables in parabolic coordinate. We have shown that the sep- 
aration procedure reduce to one-parametric (parabolic separation constant) Schrodinger type 
differential equation (30), living in the complex plane. It have been proven that the require- 
ment of convergence for solutions of eq. (30) at the singular points // — ±oo and fi = ioo 
lead only to polynomial solutions (48) with restriction for energy spectrum E in form (22) and 
at the fixed energy (or quantum number n) give the spectrum of separation constant as root 
of nth-degree algebraic equation. In difference of the previous case, the power expansions of 
polynomial solutions satisfy three-term recurrence relations and cannot be rewritting in explicit 
form and admit only numerical description. For this reason no sence to term the equation (30) 
as exactly-solvable. 

At the other hand side, the substitution of the formula of energy spectrum in the eq. (30) 
bring us to equation 



which on the real axis completely coincides for ki = 4/3a;^ and 1 ± A;2 = 25, with the one- 
dimensional spectral problem (10), and term as quasi-exactly solvable problem. Now it is 
easy to understand the birth mechanism of quasi-exactly solvable systems. The requirement of 
convergence just in real space (which is possible to determin following the [37] as the dimensional 
reduction) in singular points fi = ±oo performe together with the polynomial solutions in 
form (47) (with restriction on real or imagenery axis), which is already non complete, also non 
polynomial part. We also can shed a hght on the mistery of zeros of polynomial P„(a;^). Indeed, 
the substitution of the wave function (11) in Schrodinger equation with potential (10) lead to 
the differential equation for polynomial P„(a;^) in the same form as equation (74) (in variable 
x"^ = z), but with difference that the physical region of eq. (74) is whole real axis z G {—oo, oo), 
and therefore all zeros (for positive and negative x"^) of Pnix"^) corresponds to the zeros of 
two-dimensional eigenfunction of singular anisotropic oscillator in parabolic coordinates. 

We also shown the situation have repeated in the case of second potential (82). We have 
determined that the separation of variables in two-dimensional elliptic coordinates lead to 
Schrodinger type equation (100) in complex plane and requirement of convergence at the point 
C = 0, 27r and ( = ioo deduce to polynomial solutions and difine the energy spectrum (88). So, 
it generate trigonometric and hyperbolic quasi-exactly solvable systems (see potentials 5 and 8 
in [60]) in form of 




+ {-r + Oi{2n + 2±ki±k2) \ cosh^ v 




— cosh^ V — 
4 



i.2 _ 1 i!,2 _ 1 
^1 4 I ^^2 4 



+ A X = 0, 



sinh^ V cosh^ v 



d/i'' 




+ a{2n + 2±ki±k2)] cos^ /i 



" 4 , 

— COS II -\- 




+ A y = o. 
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where a = Thus we have proven that in the base of quasi-exactly solvabihty phenomena 

stay the dimensional reduction of superintcgrable systems on the one-dimension ^. 

This analogy prompt us to use the term of quasi-exactly solvability for the equations of 
type (30) or (100), define in the complex plane and which are not exactly-solvable but also 
admit complete set of polynomial solution. Thus we suggest to call the quantum mechanical 
systems first-order quasi-exactly solvable if the polynomial solution one-parametric differ- 
ential equation of the kind of Schrodinger equation or -dimensional equation after separation 
of variables is defined through three-term recurrence relations and the discrete eigenvalues could 
be calculated only numerically as the solutions of algebraic equations. Accordingly this defi- 
nition systems (30) and (100) are first order quasi-exactly solvable. What about quasi-exactly 
solvable systems in his old sense, it more convenient to call as partially quasi-exactly solvable, 
because only part of solution is polynomial. 

It is easy to generelize our definition for A^-dimensional quasi-exactly solvable systems. Let 
us consider the Helmholtz equation on A'"-dimensional sphere Sn- u\-\- U2 + ...u%j^i = 

^ AiB* + E^^0, (154) 
where the Laplace-Beltrami operator is 

The separation of variables in A^-dimensional ellipsoidal coordinates [82, 83, 84] 

^i-^rmrf^-^ i = 1,2, ...iv + 1 (156) 



nf=i(p,-«^) 
nf^V(a,-a,) 



with 



oi < pi < 02 < P2 < 03 < < Pn < oat+i, (157) 

lead to system of A'" - equivalent differential equations each in the region ai < pi < aj+i. They 
have the following form 

4 ^/m^i^) ^ ^m\p-a,) ^ + (^Ep^-^ + J: A,p^-^- j ^ = (158) 

and together with the energy E contain the {N — 1) separation constants Aj {i = 1, 2, ...A^ — 1) 
as parameter. This is the case of complete nonseparation of separation constants. Equation 
(158) is the generalized Lame' equation and falls into a class of equations of the Fuchsian type 
with N + 2 singularities in the points: {ai, a2, Oiv, oat+i, oo} ^. 

''Really wc can express our observation in the form of following hypothesis: all quantum-mechanical problems 
which known for today as one- dimensional quasi-exactly solvable systems could be determine via separation of 
variables in N-dimensional Schrodinger equation for superintcgrable systems. 

^"(01,02, ....ajv, ajv+i) are elementary singularities with indices (0,1/2 ) and a point at infinity is regular. 
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Consider N = 2 when eq. (158) reduce to the famous Lam'e equation. It is possible to 
prove that condition of finiteness for eigenf unctions ip in the whole region for variable p including 
the singular point (ai, 02, 03) perform only the polynomial solutions (Lam'e polynomials) which 
provide the quatization of energy in form E — i{i+l)/R'^, where i — 0, 1, 2... and simultaneously 
determine the spectrum of separation constant Xi{R) as a solution of algebraic equation. It is 
well known that Lam'e polynomials can be represented by the series expansion aroTind one of the 
singularities with the three-term recurrence relation for corresponding coefficients. Accordingly 
of our definition the equation (158) at = 2 in interval (01,03) give as the example of first 
order quasi-exactly solvable system. It is also obviously that reduction into intervals (01,02) 
or (02, 03) for the fixed values of £^ = i{i + 1)/R^ admit both polynomial and non-polynomial 
solutions. The system become partially quasi-exactly solvable. 

Let us now consider the next dimension = 3 which is more complicated. For A^ = 3 equa- 
tion (158) already contain together energy parameter E also two separation constants Ai and 
A2. As in previous case it is also possible to prove that the requirement of finiteness for wave 
function simultaneously in three intervals (oi, 02), (02, 03) and (03, 04) lead to polynomial solu- 
tions and provide quantization of energy according the formula E — J{J + 2)/R^. Accordingly 
to [67] the generalized Lam'e polynomial is derived as expansion around one of the singularities 
Oj (i = 1,2,3,4) of the equation (158) where the coefficients of this expansion obey the four- 
term recurrence relations. The obtained from this recurrence relations system of homogeneous 
algebraic equations is overcomplete since the number of equations is larger than the number of 
unknowns (expansions coefficients), and the corresponding matrix is rectangular. Concerning a 
homogeneous systems of this type, it is known that a necessary and sufficient condition for the 
existence of a nontrivial solution is equality to zero of all determinants, however, as it is proved 
in [67], for our system it is sufficient that only two determinants, resulting from this system 
by eliminating the last and the next to last rows, be equal to zero. Therefore we are coming 
to the system of two algebraic equations which give a quatization of two separation constant 
Ai(i?) and \2{R)- So, we term the equation (158) for A^ = 3 in the whole interval of variable 
p e (01,04) as second order quasi-exactly solvable. 

Now we can determine the A^th order quasi-exactly solvable systems following way: the 
quantum mechanics system is A"th - order quasi-exactly solvable if the polynomial solution 
of the A^ parametric differential equation of the Schrodinger equation type (or A^ -|- 1 - di- 
mensional Schrodinger equation after separation of variables) is defined through A^ + 1-term 
recurrent relations and the spectrum of eigenvalues of separation constant could be calculated 
only numerically as the solutions of algebraic (or systems of algebraic) equations. 

Thus in this work we classify four kind of solvabihty: exactly-solvable, quasi-exactly 
solvable, pcirtially quasi-exactly solvable, and non-exactly solvable systems. 

The interesting question which we not consider here is following: is it possible to classify 
all polynomials which coeeficients in power series expansions can be determine through three- 
term or more higher order recurrence relations as it known for classical polynomials as Jacoby, 
Gegenbauer, Legandre, Lagerre and Hermite polynomials? 
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